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ABSTRACT 


By using irregular sampling, an image may be digitized 
by concentrating sampling in areas of interest and reducing 
sampling in other areas. With this approach, storage 
required for the digital image is more efficiently used and 
is usually reduced. 

This thesis explores an interactive method for 
reconstructing a compressed _ § image. An interactive method 
allows human intuition to be combined with the speed and 
versatility of a computer. The user is able to guide the 
reconstruction through difficult cases of loop connectivity, 
branch linkup, and triangulation. MOSAIC, an element of 
MOVIE.BYU, is used for this reconstruction. 


This study shows that an interactive method has both 
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I. INTRODUCTION 


A. BACKGROUND 

Image processing can be described as the alteration and 
analysis of a picture for. such purposes as enhancement and 
recognition. Although improvements in processing methods 
for transmitted digital pictures have continued over the 
last four decades, Lt was not until the advent of 
large-scale digital computers that many of these methods 
became practical. 

Digital image processing techniques are playing a key 
role ina variety of problems. Any area that has a need for 
methods capable of enhancing pictorial information for human 
interpretation and analysis can benefit from image 
processing. The space program, an early innovator in image 
processing, enhanced and restored images from such programs 
as the Surveyor missions to the moon, the Mariner missions 
to Mars, and the Apollo manned flights to the moon. The 
improvement of the processed pictures over the unprocessed 
pictures was truly amazing [Ref. Tele In the medical field, 
physicians are assisted by computer procedures that enhance 
the contrast or code the intensity levels into color for 
easier interpretation of xrays and other biomedical images. 


Other successful applications of image processing concepts 
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can be found 1n archeology, astronomy, biology, law 
enforcement, and defense applications. 

A basic, general-purpose digital image processing system 
> snown in Figure 1.1 {Ref. 2]. The operation of the 
system may be divided into three principle categories: 

Digitization 

Processing 


Displav 


B. DIGITAL IMAGES 

A digital image is an image which has been discretized 
both in spatial coordinates and in brightness. The spatial 
resolution is the number of pixels into which an image is 
divided, indicating the precision and accuracy horizontally 
and vertically. The brightness or gray level may be viewed 
as a third axis. 

In digital image processing systems one usually deals 
with arrays of numbers obtained by spatially sampling points 
of an optical image. After processing, another’ array of 
numbers is produced, and these numbers are _ then used to 
reconstruct a continuous image for viewing. A typical 
representation of a digital image 1S amatrix. The row and 
column indices identify a point in the image and the 
corresponding matrix element value identifies the gray level 


at that point. 


ea 
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The number of gray levels and dimensions’) of the matrix 
may vary from application to application however there are 
advantages to using a number of gray levels which are 
integer powers of two and to selecting a Square array. AS a 
reference, a typical monochrome TV image isa 512 x 512 
array with 128 gray levels.. A minimum system for general 
image processing work should be able to display 256 x 256 


pixels with 64 gray levels. 


on IMAGE PROCESSING TECHNIQUES 

Techniques for image processing may be divided into four 
principle categories: 

Image Digitization 

Image Enhancement and Restoration 

Image Coding 

Image Segmentation and Representation 

Image digitization deals with converting continuous 
brightness and spatial coordinates into discrete components. 
Image enhancement and restoration concerns the improvement 
of a given image for human or machine perception. Image 
coding is used to reduce the number of bits ina digital 
image. Image segmentation and representation deals with the 
decomposition of an image into a set of simplier parts and 
the organization Simenese, palrtsamin a meaningrul cescriptive 


manner. [Ref. 2] 
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The exact approach to the above depends on which of two 
broad categories one chooses to work in. Processing 


techniques in the first category are based on modifying the 


Fourier transform of an image. The basis for this is the 
linear system convolution theorem. The spatial domain 
refers to the image plane itself, and approaches in this 


category are based on direct manipulation of the pixels in 


the image plane. 


D. PROBLEM DESCRIPTION 

Digital representations of images usually require a 
large number of bits. The logical solution is to consider 
techniques that can represent an image with fewer bits. 


When the transmission or storage of a signal requires 


excessive channel or storage capacity, the requirement can 
be reduced by more efficient coding. This is referred to as 
source or image encoding. Applications of image encoding 


fall into one of three categories: 

Image Data Compression 

Image Transmission 

Peature LxeXCraclCion 

Data compression applications are motivated by the need 
to reduce storage requirements. In image transmission 
applications interest lies in techniques which achieve 


maximum reduction in the quantity of data to be transmitted. 
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Feature extraction applications are used primarily for 
pattern recognition by computers. 

The resulting problem from the use of the above 
techniques is how to obtain an image that is acceptable for 
visual or machine analysis. Alternatively, the problem may 
be viewed as how to improve an image that has been modified 
by one of the above applictions into an acceptable image. 


This thesis will concentrate on the reconstruction of 


images that have been compressed. During the reconstruction 
phase, existing numerical algorithms are not always 
sufficient. Loop and ‘contour connectivity often causes 


problems where there are _ several loops on adjacent levels. 
This causes features on the reconstructed image to be 
distorted with sometimes unacceptable results. An 
interactive method allows human intuition to be combined 
with the speed and versatility of a computer. This should 
Besult in a better reconstructed image. A uSer can interact 
with the program to guide it through difficult cases of loop 
connectivity, branch linkups, and triangulation as well as 
to form patches. 

Programs written by CC. T. Miranda will be used _ to 
compress digital images. The specific interactive system 
used here is Movie.BYU which is distributed by Brigham Young 
University. Movie.BYU is a general purpose computer 


graphics software system. 


ic> 


The purpose of this’ study is to investigate the 
interactive procedure for image reconstruction. Answers to 
the following questions will be sought. Can this procedure 
achieve better results than other types of procedures? What 
problems are encountered when using an interactive method? 
What are the nature and characteristics of problems that are 
encountered ? What skill level is required of the 
interactive user? What hardware demands are needed to 
implement such a system? And, finally, can this procedure be 


used in practical situations? 


Iki: 


ime IMAGE COMPRESSION 


A. INTRODUCTION 

Digital representation of images usually requires a very 
large amount of data. When the number of images that must 
be retained increases, the storage requirement increases. 
An example is the storage of X-ray pictures’ in hospitals. 
Storage requirements are immense and retrieval of specific 
images is difficult. 

It is important to consider techniques’ for representing 
an image, or the information contained in the image, with 
fewer bits. Image data compression applications’ are 
motivated by this need to reduce storage requirements. An 
additional benefit gained from data compression iS a 
reduction in the quantity of data to be transmitted. With 
less data, image transmission time and the required 


bandwidth are reduced. (Ref. 3] 


B. IMAGE SURFACE APPROXIMATION WITH IRREGULAR SAMPLES 

For most images the objects of interest are localized. 
Most of the image consists of a background. This background 
has little or no variation of gray levels. As a result of 
this, using spatial sampling with regular spacing between 
the samples wastes storage space for areas with little gray 


level variation and little interest to the viewer. 


alka 


If a method could be devised that allows sampling to be 
concentrated at areas of interest and to be reduced in other 
areas, the storage required for an image could be more 
efficiently used and probably reduced. ££ £4[he preblenee.. 
to best select these irregularly spaced samples without 
causing any Silent team loss of the information in the 


image. 


Cs I B=SPLINE EUNGCATONS 


Two general approaches exist for image coding. One 
utilizes spatial domain techniques, the other operares aaa 
the transform domain of images. The approach used 


throughout this thesis 1S a spatial domain technique and is 
based on the use of B-splines. B-splines are used _ to 
approximate contours using variable knots. 

A formal definition of a spline function is as follows 
[Ref. 4]: a function S(x), defined in the Panve a Glee 
is called a spline function of degree k With Knots ce 
1=1,2,...,n if the following conditions are satisfied; 

Cr) In each interval [t,..,t, 1, W132 aw 
Cae i491 - Ds Ga 23 Se is given by a polynomial of 


degree k or less. 


(2) S(x) and its derivatives of order 1,2,...,kK-1 are 
continuous everywhere in the interval fa,b]. 


The order of the spline function is defined as 
m= k + 1 C2) 
where m is the crder of spline and =k is the degree of 


the spline. 


kes) 


The piecewise polynomial spline defined above can be 
cheracterized as a smooth function that fits at the break 
points satisfying some ending point requirements. These 
constraints have the following drawbacks [(Ref. 5]: 

OL) If we are processing experimental data points, most 
of the time we are not interested in fitting a curve 
exactly through this set of points. The reason is 
that there is an intrinsic error associated with 
these measurements. {t is preferable to fit a curve 
that presents a tradeoff between closeness of fit and 
smoothness. 

(2) If the number of data points is large the solving of 
the tridiagonal system associated with the spline 
function is time consuming and equires a lot of 
storage. 

(3) The lack of the ability to manipulate the curve 
fitting 1S a poor characteristic for a curve fitting 
algorithm. 

B-spline functions which are based on linear space 
theory overcome the above constraints and are well suited 
for use in image processing. Fundamentals of B-spline 
functions are discussed by De Boor [Ref. 6]. Andrew [Ref. 
7] discusses some early and important applications of 
B-splines for image processing. 


Some important properties of B-spline functions in the 


area of image processing include the following [Ref. 5]: 


(1) Local Basis: Only k+1 B-splines have non zero value 
at any particular interval tw, tw., J. This 
property suits very well the nonglobal 


characteristics of the pixels in an image plane. 


W2) Non-Negative Basis: The B-spline is non-negative 
which is an asset to image processing since the 
Signals represented are usually light intensities and 
are always positive quantities. 


Ire, 


(3) Differentiation and Integration: The B-spline can be 
evaluated from a recursive relation which allows the 
numeric differentiation and the integration to be 
carried out efficiently by just sudDtracting or adam 
the coefficients respectively. 


D. APPLICATION OF B-SPLINE FUNCTIONS IN NOT SELECT. 

Miranda (Ref. 5] uses the B-spline basis as “fundamental 
entities for a data compression model." Miranda’s method 
can be broken up into two steps. 

First the image contours are generated from the data 
over a rectangular and uniform grid representing the image. 
The information contained in the image signal and the number 
of contour levels or layers chosen determine the _ total 
number of contour lines. 

Next, from the image contours generated above, knots are 
selected so that the image 1s represented as data over a set 
of uneven’ knots. It contains less elements than the 
original data set, but at the same time all important 
features of the original signal are maintained. 

In order to extract important knots’ from contours, the 
contours must be stored in such a manner that 1s convenient 
for the knot selection algorithm. 

Miranda developed a series of programs that take as 
inputs an image represented as a uniform grid and user 
controlled parameters that indicate the number of levels and 


a threshold that eliminates contours with less points than 
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this threshold. The threshold level has a smoothing effect 
on noisy images. The final output of this’ series of 
programs are three files. 

File 1 is a contour record. DEmecontaumns the (X,Y) pairs 
that make up the coordinates of each ee The contours 
are concatenated together with the first two entries of each 
contour indicating the number of points on the contour and 
whether the contour is open (-1) or closed (-2). =An end 
marker (-3, -3) is used to indicate the end of data. Figure 
2.1i1is an example of a contour record. 

File 2 is a character matrix with the same dimension as 
the given image. It contains the projections of all 
contours of the image surface. Each layer is coded by a 
character starting from character A (lowest contour level) 
through character Z (highest contour level). Figure 2.2 is 
an example of a character matrix. 

ane 3 is a layer record indicates the gray ievel of 
each layer and how many contours there are at each gray 
level. Figure 2.3 is an example of a layer record. 

After the contour information has been’ stored ina 
convenient format it is time to fit a smooth spline function 
to each contour. The goal is to define a measurement for 
the smoothness of the approximating splines and another 
measurement for their closeness of fitting so that these two 
properties, usually contradictory, can be controlled by a 


single parameter. 
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Figure 2.1 Contour Record Format 
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Miranda develops and implements an algorithm that allows 
the user to automate the determination of the knots. Since 
the number of selected knots is’ much lower than the initial 
contour points, high compression rates can be achieved. 

The following example gBives an overview of this 
algorithm. For further -details see Chapter 3 of Miranda’s 
mwoerk (Ref. 5S]. 

Figure 2.4 is a superimposition of the contour record, 
Figure 2.1, and the character matrix, Figure 3.2. The 
character matrix 1s composed of fully connected contours 
represented in Vettes muhwke, the cComtounrs FreCOrd consists of 
the contours represented by unevenly spaced knots marked 
with a circle. The contour density measurement is 
accomplished by moving the isometric cross’) along’ each 
contour so that 1s is centered on every contour record point 
Merreled character in Figure 2.4). Information about the 
distribution of contours in that specific neighborhood can 
be validated by checking in four directions spaced 90 
degrees apart. The length of each arm of the isometric 
cross is set by a parameter which 1s input by the user. A 
value of five for this parameter was found by Miranda to be 
a good choice. 

The other user supplied parameter to this algorithm is 
the smoothing factor. This parameter controls the smoothing 
and fitting of the contours. A small value gives a function 


that fits closer but is less smooth due to the large number 
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of selected knots. A large value gives a smoother function 
that does not fit as close. The final choice ona value for 
the smoothing factor depends on the problem at hand and the 
desired results. 

Figures 2.5, 2.6, and 2.7 Be onsenee the effects of 
this algorithm on an image. Figure 2.5 is the original 
image that is to be processed using the programs written by 
Miranda. Using a compression ratio of 30, Figure 2.6 is the 
resulting reconstructed image. Figure Ee | is the 
reconstructed image after applying a compression ratio of 


40. 
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Sie IMAGE RECONSTRUCTION 


A. BACKGROUND 
The true test of any image coding algorithm must _ be 
based on a comparison . of the original image with a 


reconstructed image. The criteria for this comparison may 


vary from application to application. Depending on the 
fcended use of a reconstructed image there may be 
acceptable differences between the original and 


reconstructed images. 
There is no single optimum method for image 
reconstruction. There is a "best" restoring method relative 


to the type of a priori information regarding the object and 


noise one might have. Other, more peripheral factors, such 
as the amount of data to be processed, and permissable 
computer cost, also must enter into the choice. These 


factors must often be balanced against an "optimum" choice 
based purely on accuracy. [Ref. 1] 

A major issue in image reconstruction is’ tne intended 
“evaluator” of the reconstructed image. A person may pass 
favorable judgement on a reconstructed image while a machine 
may reject the reconstructed image. A person basically 
makes a qualitative judgement between the original and 
reconstructed images. On the other hand a computer usually 


makes a quantitative evaluation of the two images. This 


Sel 


quantlvratlive evaluation is based on mathematical 


relationships. 


B. QUALITATIVE FIDELITY CRITTER. 

When the output images are to be viewed by people who 
usually use a subjective fidelity criteria CoOrrespondingwre 
how good the images sneer to human observers, it is best 
to use qualitative Criteriae 

Since digital images are displayed as a discrete set of 
brightness points, the ability of the eye to discriminate 
between different brightness levels iS an important 
consideration in presenting image processing results. The 
human visual system can adapt to ae wide range of light 
intensities. This range is on the order of 10'° from the 
scoptic threshold to the glare limit. 

A key point however is that the visual system Can not 
Operate over this entire range Simultaneously. rt 
accomplishes this large variation by sliding the center of 
its sensitivity range, a phenomenon known as brightness 
adaptation. The human eye 1s only able to differentiate 
between about 30 gray levels and about 120 different colors. 
[Kets 

Another feature of the human eye 1s its tendency to 
filter out any abrupt changes in brightness levels. This 


can be compared to a smoothing effect. 
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There are no mathematical relationships that can 
indicate the qualitative fidelity of a reconstructed image. 
It is entirely a subjective matter. A commonly used scale 
Home CepOorting the qualitative fidelity of an image was 
developed by Panel 6 of the Television Study Organization 
(Ref. 8]. This scale ranges from excellent to unusable with 
the following definitions: 


mee excellent: The image is of extremely high quality, 
as good as you could desire. 


(2) Fine: The image is of high quality erov 1 dine 
enjoyable viewing. Interference is not 
objectionable. 





(3) Marginal: The image is of acceptable quay 
Interference 1s not objectionable. 


(4) Inferior: The image 1s of poor quality but you 
could watch it. Objectionable interference ls 
definitely present. 

(5) Unusable: The image 1s so. bad that you can not 
Watery it. 

SeeeeCUANTITATIVE FIDELITY CRITERIA 

Some image transmission systems can tolerate errors in 
the reconstructed image. When this 1s the case, a fidelity 
criteria can be used aS a measure of system quality. The 
basic measurement of the error between an input pixel and 
the corresponding output pixel is the difference in gray 
levels between these two pixels or expressed mathematically 

etx,y) = gtx,y) - flx,y) (Sat) 

where g(x,y) is the output image signal and (x,y) is the 


spatial coordinates. [Ref. 2] 


oS 


One can apply this error  tfunetion, 8a. in several 
ways to come up with a measurement of quantitative fidelity 
criteria. Some examples of this are the root-mean-square 
(rms) error between the input’ image and output image, and 
the rms signal-to-noise ratio of the output image. The 
root-mean-square error 1s defined as the square root of the 
squared error averaged over the entire image. The 
signal-to-noise ratio can be defined as the square root of 
the peak value of g(x,y) squared and divided by the 
root-mean-Square error. Gonzales [Ref. 2:pp. 229-2311] has a 


complete development of the specific equations. 


D. SURFACE RECONSTRUCTION 

Now that some guidelines for comparing a reconstruction 
image with the original image have been laid, a method must 
be found to reconstruct the set of irregular samples given 
DY “Xe. Var; ee ae le eee ae ee Is 1s necessary to f imeia 
smooth function F(x,y) so that F(x, , Yud = ff. , K=d5 2 eee 
Many interpolation methods exist that are based on either a 
global interpolation method or a local interpolation method. 

To simplify the reconstruction process and Keepeine 
computational time required reasonable, local interpolation 
methods will be used. 

In local interpolation methods, F(x,y) depends only on a 
limited neighborhood of § f, ’s. The justification for using 


these methods is that the correlation between two pixels in 
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an image usually decreases rapidly with the distance between 
them. The influence a sample exerts on another sample is 
inversely proportional to the distance away from it. Two 
eommoniy used local interpolation methods are the triangle 
element method and the local thin plate spline method. 

The method to be used in this paper is based on the 
triangle element method. In the triangle element method, 
the first step involves building ae triangular network 
covering the convex hull of the samples. The (x. ,y.) points 
become the vertices of the triangular cells. The second 
step involves the Evaluation ob an arbitrary point (x,y) in 
the output array. This is done by linearly interpolating 
Maen pOolmne from the values at the vertices of the triangle to 
which the point belongs. A step by step procedure to 
construct a uniform grid over ae set of scattered points 
(x,y) 1s given by Lawson in a comprehensive paper [Ref. 9]. 
Figure 3.1 1S an example of a triangulation grid. 

Unfortunately any partition of a convex hull WE, 
triangular cells is not unique. some triangularizations 
wee be Superior to others in the sense that less long thin 
triangles are generated and a better reconstruction of the 
Original image results. 

The algorithm discussed by Lawson is’ well suited for 
automatically generating these triangles by computational 


methods. A brief summary of the method follows. 
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Mealro (Uia(etaw a aed Triangulation Grae 


Using the three closest pairs of points, elie iets St 
triangle is constructed. A point is added ina sequential 
manner in the ascending order of the distance from the 
midpoint of the existing edge visible to it. The added 
triangle and tne neighboring triangle together form a 
quadrilateral. The algorithm next tests whether the 
alternative dissection Can maximize the minimum interior 
angles of the resuitant two triangles. Whenever necessary, 
an exchange of diagonals in the quadrilateral is conducted. 

Unfortunately automatic numerical algorithms are not 
always sufficient in the reconstruction = phase. Several 
hagins OM adjacent ~contour levels w~often causes a problem 
where triangulation does not follow the contour lines. This 
causes features on the reconstructed image to be distorted 
with sometimes unacceptable jagged contours. By being able 
to interactively guide the triangulation process”) through 
high gradient areas, the reconstructed image will hopefully 
show an improvement over the image reconstructed by a method 
Similar to Lawson’s. 

The remainder of this paper discusses an interactive 
method for triangulation of a compressed image and compares 


results with the triangulate method. 
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4. MOV EERE YU 


Ae VER V4 EW 

In order to interactively reconstruct an lmage one wngee 
have a set of software tools. This is not a simple matter. 
The complexity required of such a program is large. This is 
especially true when dealing with images that include both 
man-made objects and natural objects. Man-made objects 


usually have well-defined shapes that are easily discernable 


when they are viewed in contour lines. The opposite is true 
for natural objects. These contour lines do not always give 
a true indication of the features they represent. In the 


presence of noise their irregular shapes usually cause 
problems in reconstructing the image for display. Any 
program tools used to interactively construct an image must 
have a way to allow the user to deal with this kind of bad 


contours. 


Be INTRODUCTIONS TO. MOVIE. BU 

The particular set of programs selected for use in the 
image reconstruction process is MOVIE.BYU Version (5). 
MOVIE.BYU is a general purpose computer graphics software 
system distributed by Brigham Young University. The basis 


for the remaining discussion on MOVIE.BYU is the MOVIE.BYU 
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Training Text (Ref. 10] and the MOVIE.BYU User’s Manual 
PRefo- it). 

The elements of the MOVIE system are FORTRAN programs 
for the display and manipulation of . data representing 
mathematical models. The geometry of these mathematical 
models may be described in terms of polygonal elements, 
polyhedral solid elements, or contour line definitions. 

All input to the MOVIE system is~ read as alphanumeric 
characters with each character interpreted individually. 
Tinis™circumvents all of the restrictions and conventions of 
FORTRAN input formats and makes the code virtually "bomb 
preoer.' = Phe key words, integers, and real numbers are then 
reconstructed in software and the commands are interpreted 
from the resulting list of key words and numbers. The 
current edition of MOVIE.BYU has also been enhanced with 
Capabilities to generate and access command files. This 
Simplifies repetitive input of the same series of commands. 

The hardware requirement for MOVIE.BYU is a time sharing 
digital computer with a word length of at least $32 bits. 
Software interfaces for most major graphic terminals are 
available including a special software interface for the 
Tektronix 4027 which Significantly extends the color 
possibilities by the use of patterns which’ simulate 


additional intensities for each color gun. 


ene, 


C.- ELEMENTS OF MOVE Snal 

The MOVIE system is composed of six elements - DISPLAY, 
UTILITY, SECTION, TITLE, MOSAIC, and=3GeMreee "The six 
programs in MOVIE work in harmony to prov.ide displays of the 
data in line drawing or continuous tone image format, to 
clip and cap three dimensional systems to expose internal 
surfaces, to modify geometry, displacement, and/or scalar 
function files by way of correction, appendage, or symmetry 
operations, to generate new models or representations, to 
convert complex contour line definitions into polygonal 
element mosaics, or to view multiple models simultaneously.” 
(Ref. 1i:pp. 1.1]. A brief description of each element 
follows. 

DISPLAY is the “heart” of the system. It 1s ‘ah 
interactive program for the display and animation of any 
model composed of polygons. DISPLAY allows the user to 
manipulate the model (rotate, translate, ete. specify 
colors for the background and the different element parts, 
and select the display device. 

UTILITY is a data generation and editing program which 
allows the user to produce and/or edit models of two or 
three dimensional polygonal systems. The FORTRAN data files 
created are ina format which 1s compatible witn the other 
programs in the MOVIE system. 

SECTION 1s a program used to modify solid data 


representations so that they are compatible with DISPLAY. 
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TITLE is a program to generate two and three dimensional 
characters whose data format is the same as that used Dy the 
other programs. The user enters a line of text which is 
converted into characters composed of polygons according to 
specifications given by the user. 

MOSAIC is an interactive program that converts complex 
surfaces defined by contour lines into mosaics of triangular 
and quadrilateral elements. 

COMPOSE is a program that allows’ the retrieval of data 
penerated by DISPLAY (with the RECOrd option) to compose 
drawings with multiple images and text. 

Of the six modules listed above only two are of interest 
in this work. MOSAIC is the program that allows a user to 
interactively control the triangulation of a muliple level 
contour representation of an image. DISPLAY is of interest 
in that it allows the user to get a three dimensional 
perspective of the connectivity of the layers after the 
triangulation 1S completed by MOSAIC. The user can then go 


back to MOSAIC to resolve any inconsistencies’ that showed 


oD. 


ee MOSAIC 
ihe Paeroauctl ion 
MOSAIC implements) an algorithm for processing 
complex contour arrangements into polygonal element mosaics 


which are suitable for line drawing and continuous” tone 
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displays. The program does this by mapping adjacent 
contours onto the same unit square and, subject to ordering 
limitations, connecting nodes of one contour to their 
neighbors in the other contour so that the total length of 
the connecting lines iS-— minimized. A user 1S able to 
interact and resolve highly ambiguous situations. Other key 
features of MOSAIC are node thinning and selective reduction 
of triangular pairs to quadmadavera lice 

Figure 4.1 illustrates a surface which has_ been 
reconstructed using MOSAIC. Figures 4.1a and 4.tc show the 
original contour loops prior to processing, Figures 4.1b and 
4.1d show the triangulation performed by MOSAIC. 

2. Triangulation Algorithm in MOSAIC 

The triangulation algorithm implemented by MOSAIC 
assumes that the optimum arrangement of triangles produces 
the least sum total length of all connecting diagonals 
between levels. MOSAIC also requires that the loops to be 
triangulated are numbered in the same rotational directaGgnm 
and that the first nodes of each loop are closed. This is 
something, however, that the user need not worry about since 
MOSAIC automatically performs these preparations. 

This triangulation algorithm is illustrated in 
Figures 4.2 and 4.3. Figure 4.2 depicts two loops that are 
to be triangulated. The graph shown in Figure 4.3a 1S a 
representation of these loops. Each intersection of lines 


can be thought of as defining a diagonal, stich that Gia. 
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depicts the diagonal connecting node i on the bottom coneoun 
to node j on the top. Associated with the graph of Figure 
4.5a are the two matrices shown in Figure 4.3b and Figure 
4.3c. Both the total distance matrix and the path memory 
matrix coefficients have a direct correspondence with the 
points on the graph. 

The total distance matrix is the sum of the length 
of the diagonals included in the shortest path. The path 
memory matrix is used to indicate the path that was used to 
attain the least total length in lacing from diagonal P(1,1) 
bOur Ci pe M(i,j3) equals one if the path including diagonal 
ib-jt also includes the diagonal Cie Web 7 te If the path 
includes diagonal ib-(j-1)t instead of (1-1)b-jt, M(1,)) 
equals zero. 

For best results when using this’) algorithm  toop 
pairs should be of similar size and shape. The loop pairs 
should not be offset excessively. For this reason, MOSAIC 
automatically maps both loops onto the same unit square 
before triangulating. 

oO. Loop Connectivity and Branching 

When more than one loop exists on adjacent levels 
potential problems may occur. To alleviate this problem 
MOSAIC allows the user to interact and specity which Woegee 
on the bottom layer connect to which loops on the top Tavern 
As a default MOSAIC uses an overlap test to determine which 


loops to connect. See Figure 4.4 for an example. This is 
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good for many cases but there are cases of COnneCtIvitye ioe 
Overlap will not resolve correctly. After the user has 
forced any desired loop connectivity he may return the 
control of the process back to the algorithm. 

The other problem that sometimes occurs when there 
are more than one loop ona level is branching. The method 
used by MOSAIC treats all branches on a level as one 
continuous loop. It does this by locating the closest nodes 
between branches and then renumbering the branching loops so 
that they collectively appear to be one large loop that has 
been squeezed eee tate over segments where the  1linkups 
OCCUR. It then triangulates as before. If the user is 
unhappy with this branching he may force branching between 
given nodes of the loops. Figure 4.5 illustrates | Gime 
concept. For the loops shown in Figure 4.5a, a flat patch 
between the two co-planar loops would be a better solution 
than the automatic routine. To accomplish this, the user 
must input the node numbers associated with nodes A, E, B, 
and D. MOSAIC will then automatically form the patch by 
lacing the two loops together from nodes A to E and B to D. 
The remaining portions of the top loops are then renumbered 
and triangulation carried in a normal manner. 

4. Other User Options 

A large data base of contour information sometimes 

becomes difficult to use. If this is the case, it may be 


desirable to eliminate unessential nodes. MOSAIC has four 
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different ways in which it can thin out a data base. Peis 
able to eliminate nodes that are too close together. It ‘cam 
combine line segments between which the change of direction 
is minimal. [t can skip entire contour levels as the data 
is read in. And finally, it can join pairs of adjacent 
triangles to form a quadrilateral if the angle by which they 
are out of plane is sufficiently small. 

All of these options require input by the user 
regarding to the acceptable distance between nodes, amount 
of direction change between elements, 2 level spacing, and 
quadrilateral warp. 

5. Command Structure 


MOSAIC uses a multi-level prompt command structure 


fOr User Inveract ion: There are a total of four levels ane 
one pseudo level. At each level the user is either asked a 
specific question or given a prompt (")>" for level 1) )e8 
for level 2.,°7>>> * for levek —4. “Oe for level 4, and 


"D2 Tor, eve l*5) 2 At any time the command HELP may be 
entered to obtain a complete listing of all availapie 
commands for the current level. "Escape" to a lower 
(numbered) level is accomplished by entering a carriage 
LeGuUrh: 

In addition to the commands unigue to each level 
there are eight global commands which may be entered at 
command levels one through four. These commands are DEVIce, 


EXIT, HELP, MAP, PART, TOTAls, WARP, and WRiTe°> foresees 
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commands and all other commands only the first four letters 
need to be entered by the user. The MOVIE.BYU User’s Manual 
(Ref. 11] contains a complete description of commands, and a 


brief description of each command is contained in Appendix 


A. Table I is a listing of all available commands. 
Seeinstallation cheuesa IC 


Using Chapter 8 of the MOVIE.BYU User’s Manual (Ref. 
113), MOSAIC was installed on a VAX 11/780 computer running 
under a VMS operating system. The version of VMS was 4.2. 
The language used was FORTRAN 4.3. 

Due to the large size of the contour data Dase to be 
used with MOSAIC, it was necessary to increase most of the 
parameters associated with MOSAIC. It was also necessary to 
increase the dimension(s) of the corresponding arrays 
associated with these parameters. 

Table II lists all parameters changed along with 
their original values and current values. Table III lists 
all arrays whose dimensions were changed along with their 
dependent parameter, original dimension, and current 


dimension. 
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TABLE I 


Commands Available In MOSAIC 


COMMAND LEVE Ee COMMAND LEVEE 


DEVIce EDIT-DELEte node 
EXLE EDIT-LOOP 

HELP EDIT-MOVE node 
MAP INSPect 

PART MANUal 

TOTAIs MANUal-AUTOmatic 
WARP MANUal-BRIDge 
WRITe MANUal-CLEAr 
CLIP MANUal -DRAW 
CLOCkwise MANUal-GUIDe 
CONVert GUIDe-AUTOmatic 
MAKE GUIDe-DENSity 
RANGe GUIDe-ERASe 

READ GUIDe-RENUmber 
THINning GUIDe-TRIAngulate 
AUTOmatic MANUal - INCLude 
CAP MANUal-NEXT 

EDIT POST 


EDIT-ADD node 
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TABLE. al 


Changes To MOSAIC Parameters 


PARAMETER | EXPLANATION ORIGINAL | CURRENT 
VALUE VALUE 


maximum number of 
loops 


maximum number 
per level 


maximum number 
per level 


MNODES maximum number 


MNPLP maximum number 
per loop 


NBNMAX maximum number loops 
displayed 


NPTMAX maximum number of element¢ 
NVMAX maximum number of nodes 


being triangulated on 
adjacent levels 





oye 


TABEE Si 


Changes to MOSAIC Arrays 


ARRAY DIMENSION SUBROUTINES USED IN 
PARAMETER 


Da hn @ a NVMAX ATRIAN, BKTRAK 

P2 NVMAX ATRIAN, BKTRAK, INIMAP 

V NVMAX ATRIAN, AUTO, BKTRAK, DFINV, 
DRAWV, FPATCH, GUIDE, INIMAP, 
LABLP, LINK, MAKV, MANUAL, MAPP 
REORDR 

TEMP NVMAX REORDR 

FLAG MLP AUTO, CONECT, DFINV, FRPATCuE 
LABLP, LINK, LNKMAN, LPLST, MAKV 
MANUAL, SELECT 

EX MLP BLOCK, CAP, CLIPIT, CONECT, Drie 
DRAWLP, EDIT, GEOMWR, LABLP, 
LNKMAN,LPLST, MAKV, MANUAL, 
PROCES, THNOUT 

P1 MLP BLOCK, CAP, CLIPIT, CONECT> Uris 
DRAWLP, EDIT, GEOMWR, LABLP, 
LNKMAN, LPLST, MAKV, MANUAL, 
PROCES, THNOUT 

LIST NBNMAX DRAWLP, EDIT, LCABLPRSErPiEst 

NIPL MNLPL ATRIAN 

LPSTK MNLPL AUTO, CONECT, FGINV, FEATGHE 
LABLP, LINK, LNKMAN, LPLST, MAKV 
MANUAL, SELECT 

O MNLPL AUTO, DFINV, CINK, ENKMAN ee 

LBN MNLPL BKTRAK, DFINV 

NPL MNLPL CLIPIT, GEOMRD, PROCES Stole 

C MNLPL CONECT, FPATCH, LPUST) s4AKve 
MANUAL, SELECT 

CP MNLPL DFINV 

D MNLPL DF INV 

Dal MNLPL DFINV 

ORDER MNLPL DFINV 

LCF LG MNLPL DFINV 

BRG MNLPL DFINV, FPATCH, LINK, LNKMAN 

NCF MNLPL MANUAL, SELECT 
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TABLE i ils Ceon: 1nued) 


Changes Made to MOSAIC Arrays 


ARRAY DIMENSION SUBROUTINES USED -[N 
PARAMETER 


BRANCH MNLPL pe Lec i 
STACK MNLPL SELECT 
EPNTI MNPLP CADE TE 
NCP MNPLP Cri PiT 


eri P MNPLP Cole it 


5 MNODES Pou AN ewan PRAKe., CAPl= CLIEPia: 


CEKWSE, CROSoeeDFINV, GEOTWR, 
UA ee otemmtaeP, MOVE, PREPLT 
THNOUT 


PS MNNPL CLKWSE, GEOMRD, MAKE, THNOUT 


ce NPTMAX ATRIAN, BKTRAK, CAP, DRAWEL, 
GEOMWR 
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V. EXPERIMENTAL RESULTS 


A. PRELIMINARY WORK 

After installing MOSAIC, the first question )tomia— 
answered was what images to use for the reconstruction 
process. The choice was quickly narrowed down to compressed 
images with either a compression ratio of 30 or 40. The 


reason for this was to have an image that had easily 


noticeable differences from the original image. It was 
finally decided to use a compression factor of 40> 50Fon 
reasons to be discussed later this was a good choice. The 


digitized image of Figure 5.1 is the basis for the remainder 
of this study. 

Next it became necessary to convert the format of the 
digitized compressed image into a format suitable for 


MOSAIC. MOSAIC requires the data to be ina single file 


that contains” the number of nodes in a loop, “thew 
coordinate of the contour level. and then the Oe. 
coordinates of the nodes in the loop. This information is 


required for each loop. A "O" at the end of the last loop 
terminates the file. The Z coordinate corresponds’ to the 
pray level of the contour. 

The information needed to construct the input file for 
MOSAIC is contained in two separate files. These two files 


are the Contour Record and the Layer Record which were 


one. 


ar, 


opty! 





Bis (Chavon oat Original Compressea Image 


fy, 


discussed in Chapter 2. To accomplish this conversion the 
program CONVBYU.FOR was written in FORTRAN. The source code 
is contained in Appendix B. CONVBYU.FOR accomplishes its 
purpose by opening and reading the two files, Contour Record 
and Layer Record, and then merging the information needed by 


MOSAIC into one file. 


B. RECONSTRUCTION USING MOSAIC 
1. Closed Contours 

With MOSAIC installed and the contour data needed by 
MOSAIC converted to the correct format, MOSAIC is executed 
and the contours are read in. Using the INSPect command all 
levels were displayed individually. The first thing noticed 
was the fact that all contours were closed, even those that 
were originally open contours. It turns out that MOSAIC 
requires all contours to be closed and, if a contour is not 
closed, MOSAIC closes the contour by connecting the first 
and last nodes. 

By closing contours ine Jenis arbitrary manner, 
contours) that were created sometimes crossed over 
themselves. This could prove confusing to the triangulation 
algorithm in MOSAIC and result in a poor reconstruction. 

After further study of the individual levels, 1 
appeared that instead of closing each loop a better choice 
might be to combine several groups of loops into one loop 


per group. 


a3 


A quick and easy wasy to Simulate the combining of a 
pair of loops and achieving the desired appearance of one 
loop in triangulation is to use the BRIDge command. 
Unfortunately due to the manner in which the loops were 
closed by MOSAIC (contour crossover) this was not a viable 
option. 

The only other option was to use the EDIT commands. 
By adding, deleting, and/or moving nodes, two contours could 
be merged intO one contour. This was an easy but time 
consuming task. It is necessary to duplicate the nodes of 
loop one on loop two and then go back and delete the nodes 
of loop one. 

The choice of which loops to combine has no set of 
rules. It comes down to a decision made by the observer. 
This decision is based on a_ study of the loops ona level 
and the visual relationship of the loops on adjacent levels. 
When contours are fairly large and have a well-defined shape 
it 1S an easy matter to determine which loops to combine. 
It is a difficult endeavor when dealing with small loops 
that have no well-defined shape. 

Elepures S.2 through 2.17 show the “original contour 
representations of the eight levels of the digitized image 
and the modified contour representations obtained through 


adding and deleting nodes as discussed above. 
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2. Loop Conmectiwaty 

With all contour levels redrawn, decisions must be 
made on which loops to connect to which loops. When left to 
proceed automatically MOSAIC uses an -overlap criteria to 
decide which loops to connect. Any lower loops that branch 
onto any upper loops will be triangulated together. 

This works fine for contours that have not suffered 
any degradaticen of Shape due to noise, transmission loss, or 
compression of the original image. When a contour’s shape 
has been degraded, it sometimes results in a loop bDranching 
onto another loop. MOSAIC will then attempt to connect 
these loops together. 

The alternative 1s to specify specifically Whiten 
loop(s) to connect to which Iloop(s). Once again this isa 
Subjective decision based on an examination of adjacent 
levels and formation of a mental three dimensional picture 
of the levels. With well-defined contours this 1s not 
Giit teu ey Other contours PFESEnt difficultiess 
Additionally, this stage of triangulation relies on accurate 
choices being made on how the contours were closed. 

Once all levels have been triangulated DISPLAY, 
another module of MOVIE.BYU, may be used to view ona 
display device a three dimensional representation of all 
levels along with the connections”) made. If this shows 
unsatisfactory results, the user may then’ go back into 


MOSAIC and respecify which loops CoNconneer 
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Although this process is Por the most part 
Subjective, they are a few guidelines to aid the user in 
deciding which loops to connect. 

Triangulating loops on two levels can be compared to 
"sewing". With this in mind one would expect’ the 
triangulation of loops to be concentrated along the edges of 
the contours. The interior of the inner loop should be 
relatively empty. When the triangulation causes lines to be 
drawn back and forth across the contours unevenly it is an 
indication that something is not right. The cause could be 
that loops were not closed correctly or that the wrong loops 
were connected. A user may not always be able to resolve 
this inconsistency. This results ina choice of what loops 
not to connect to other loops. 

Once the user is satisfied with the triangulation 
process, MOSAIC allows the user to save this data to a file. 
The file contains the (X, er Zi) coordinates of all nodes 
along with the node numbers of the vertices of all triangles 
formed. 

Figures 5.18 through 5.24 show the triangulation 
between each pair of ea acon levels that was’ finally 
settled on. 

3. Reconstruction of the Image 

The reconstruction process of the image consists of 

two steps. The first step is the triangulation process 


discussed above. The second step involves evaluating an 
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Figure 5.24 Triangulation - Loops 7-¢ 


@eottranry pOImt (xX, YY) Of a uniform (256x296) grid to 
determine which triangle it lies in. After the appropriate 
triangle is computed, a gray level is assigned to this point 
by doing a linear interpolation inside . the triangle where 
fae pDOint lies. 

A program, FRA.FOR, was written to accomplish this 
step. This program reconstructs the image as an uniform 
Pama Or pixels with the same size as the original (256x256). 
To accomplish this the data file created by MOSAIC is fed 
into FRA.FOR. The output of FRA.FOR is a uniform (256x256) 
grid representing the reconstructed image. Source code for 
FRA.FOR is contained in Appendix B. 

The COMTAL system was used co display this 
reconstructed image. Before this could be done the uniform 
(256x256) grid created by FRA.FOR must be converted into a 
COMTAL format, a uniform Col2xS12) grid. DISPoie FOR 
accomplished this by carrying out a two by two expansion. 
The source code for DISP3M.FOR is contained in Appendix B. 

4. Analysis of the Reconstructed Image 

Figure 5.25 is the original compressed image = and 
Figure 5.26 is the reconstructed image. A qualitative 
analysis of the two images reveals several points. Farset, 
in the original image the edges of the contour levels do not 
appear sharp. The jagged edges appear to be a series of 
Z’s. Comparing this to the reconstructed image, the edges 


are more well-defined with a low overshoot. 
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On the minus side there is a loss of detarieineene 
reconstructed image where the image consists of small, 
irregular shaped contours. This can be attributed to the 
difficuly in determining how to close these loops properly 
and to which loops they connect to. Also MOSAIC is unable 
to handle any loop that has less than three nodes. image 
meant that these loops were not read into MOSAIC and 
therefore had no effect on the reconstructed image. 

The final decision of whether the reconstructed 
image can be considered an improvement over’ the original 
image can be AEE? walled favorable. The major features 
of the image are definitely improved. The degradation of 
the small detail is something that the normal viewer would 
not notice. 

Next a quantitative analysis was done comparing both 
the compressed image and the reconstructed compressed image 
to the original uncompressed image. 

The definition of error used was the root mean 
square error. It was chosen since it is probably the most 
commonly used method of error measurement used in image 
processing. Background on RMS rror was discussed in 
Chapter 3. RMS error may be defined as the square root of 
the squared error averaged over the image array. The 
averaging sense of this definition resembles that Of §eie 
intrinsic human visual system. It tends to filter out high 


frequencies and introduces smoothing effects in the image. 
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The program RMSER.FOR was used to compute the RMS 
error. Source code for this program is included in Appendix 
B. 

The results of this analysis showed an RMS error of 
17.3 percent for the original compressed image and 24.8 
percent for the reconstructed image. This 1S not surprising 
Since before even starting to use MOSAIC for triangulation 
all contours with less than three nodes were lost due to 
limitations of MOSAIC. The second factor that contributes 
Meetnis error rate is the fact 
meat Other small contours were effectively ignored because 
their connectivity to other contours could not be 
determined. Still the percentage of error is relatively low 
when compared to the error rate of the original compressed 
image. 

mil 2A | Gi SU gee | contains a flow chart for the 
experimental work carried out. The next chapter will 


discuss lessons learned and recommendations for future work. 
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VI. CONCLUSIONS 


nee) 6COMMENTS 

Results obtained in this study can be summed up as both 
promising and disappointing. There is no clear outcome one 
maiyeor the other. "On the positive side the user is able to 
improve the presentation of contours that represent large 
man-made objects or well-defined terrestial objects. The 
edges of the contour levels represented in such images are 
sharp with less crossover into adjacent levels. 

On the negative side, the interactive restoration scheme 
used here results in ae significant loss of detail in the 
information contained in small, irregular shaped contours. 
Depending on the image being reconstructed or the use of 
this reconstructed image, this may or may not be of 
importance. | aa the image 1s intended for visual 
mterpretation only, the loss of detail may not even be 
noticed. However, if the image is to be subject to computer 


analysis, it may make a big difference. 


B. LESSONS LEARNED 
ite MOVIE BYU 
Although MOVIE.BYU is ae very powerful set of 
software tools, it may not be the ideal set of programs to 


use in image reconstruction. MOSAIC is better at 


sik 


manipulating contours) that are well-defined andere 
naturally. When it iS given irregular contours’ that are 
manually closed, undesirable results sometimes arise after 
the triangulation algorithm is applied. 

Some algorithms overcome this by triangulating a 
surface based on the (X, Y) coordinates of the nodes with no 
attention paid to the gray contour levels. This, hG@wevere 
Zlives rise to jagged edges cf pray levels when the 
reconstruction of an image is finished. 

Another minor problem with MOSAIC is the amount of 
work needed to install MOSAIC. As originally configured: 
MOSAIC 1s capable of only handling a esmall number of 
contours with a small number of nodes. This is “nee 
sufficient to handle the complexity of most images. To 
modify MOSAIC to handle these images one must not oni 
increase the values of the appropriate parameters but also 
search through all subroutines and increase the dimensions 
of arrays that are related to these parameters. 

A more serious problem is the fact that MOSAIC is 
not a very forgiving program, Entry of a wrong number can 
cause repercussions Shcouer all contour levels and not just 
the contour being currently worked on. Correction Of SG 
error often becomes a_ difficult, if not impossible, task. 
Unless the error can be discovered and corrected quickly, it 
is usually more time efficient to start over from the 


beginning. The user soon learns after a couple of mistakes 
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to be extremely careful when entering data or risk losing 
several hours of work. 

In its defense, MOSAIC is a very complicated program 
that offers a large degree of user interaction. Due to the 
Semelexity of the above problems there 1s no easy solution. 
When dealing with irregular shaped contours it is a very 
difficult task to offer solutions to every possibility. 

2. Hardware Requirements 
As the complexity of a contour level increases it 


becomes more difficult to distinguish individual nodes and 


loop numbers. This has a detrimental effect when it comes 
fommc losing contours manually and connecting contours 
together. 


The Tektronix 4010 proved to be no more than an 
adequate display device. This is easily confirmed by 
examining the contour levels presented in Chapter V. It is 
impossible to accurately read the loop number of every loop. 
Better reconstruction results may be expected if the display 
presentation of the contours is improved. 

5. Skill Level of the Interactive User 

The skill level of the interactive user plays a Key 
role in the reconstruction process. This is because 
interactive reconstruction methods are of a subjective 


Nature. 
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A user must be able to spot inconsistencies’ in the 
EONtCOUr Dresencataous-: Unfortunately there 1S “no sam 
encompassing list of possible errors. The user must catch 
these inconsistencies through a combination of general image 
processing knowledge and common sense. AS an example, image 
photographs are not able to show whatmienbele nme layer, they 
only show what is on top. The user should be able to notice 


when the triangulation process’ starts to violate this basic 


fact. 

With a basic understanding of image processing, the 
USer ” Can. be. -expecced We progress in his) ability to 
accurately reconstruct an imagee He will be able to spot 


prospective trouble spots through the experience he jee 
gained. Also his understanding of the effect of available 
program commands will allow him to manipulate contours 
easier. Although a new user is— able to build on the 
experience gained by others, he must still progress through 
a learning curve of his own before he can be considered 


DEOTITCliemt-. 


CC. CONCLUSIONS 

There is a need for interactive restoration techniques. 
Due to the complexity of images there will always be 
Situations where a conventional restoration algorithm falls 
short. This is especially true when the primary use of a 


restored image is for visual presentations. 
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More work must be performed in the area of 
reconstruction of small, irregular shaped contours. At the 
present time it takes a combination of a highly skilled user 
and prior Knowledge of what to expect from the contours to 
obtain an acceptable reconstruction of such a region. 

Interactive programs must be developed that’ give the 
user even more control in the manipulation of contours. A 
possible enhancement is the use of a trackball or mouse to 
control the addition, deletion, or movement of nodes. 
Another essential feature of a program is the ability to 
Beem im ON a particular area of interest. This allows the 
user to pick loop and node numbers’) of interest when an area 
1s congested with information. Both of these features were 
lacking in MOVIE.BYU. With these features it is felt that a 
better image reconstruction could be obtained. 

Additional research in the area of image combining may 
prove of interest. The word "“optimum"” as applied to 
interactive reconstruction usually implies an optimum 
response of the human visual system. When "optimum" is used 
in conjunction with automated algorithms, it usually is in 
mererence to a mathematical concept. By combining the best 
features of each image obtained with different approaches 
into one image, it should be expected that a superior 
reconstruction result will result. 

Image processing iS an area that has experienced 


vigorous growth recently. It can be expected that this 


o> 


ProOw ciiew 1) eCOnel nue: Along with this’ growth will be an 
increased need for accurate reconstructed images. With 
improved software and display devices interactive methods 
will certainly play a role. The reason is that, although 
machines play a significant part in image processing, it is 
a human observer who must ultimately interpret or approve an 


interpretation of anmunage. 
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APPENDIX A 


MOSAIC COMMANDS 


A. GLOBAL COMMANDS 


1. DEVIce allows the user to select the display device. 
2. EXIT provides a controlled exit from the program. 
3. HELP types the available commands at the Current 


point on the user’s terminal. 

4. MAP functions as a toggle switch and allows the user 
to bypass the mapping algorithm or to restart its use. 

5. PART allows the user to group elements into parts as 
they are generated. 

S. TOTAlsS causes the number of parts, nodes, elements, 
and the total surface area of the elements to be typed on 
the output device. 

7. WARP allows the user to specify the angle which 
controls’ the combining of triangular elements into 
gGuadrilateral elements. 

8. WRITe enables the MOSAIC format contour data to be 
written to disk, the DISPLAY format panel data to be written 
to disk, and the DISPLAY format contour line segment data to 


be written to disk. 


of 


B. LEVEL ONE COMMANDS 

1. CLIP acts as a toggie switch. When the switch is 
turned on the program requests the data necessary aaron 
elipping: The clipping plane algorithm can accommodate any 
number of clipping planes. 

2. CLOCkwise allows the user to specify the ordering of 
the contour loops. 

3. CONVert proceeds to level two. It has no value 
unless some data has been read. 

4a. MAKE allows the user to generate a contour data 
file. 

5. RANGe allows the user to specify which portion of 
the contour data file is to be read into MOSAIC. 

6. READ reads in the data file and prepares it for 
triangulation. If the CLIP, CLOCkwise, or THINning commands 
are to be used they must be invoked prior to the issuance of 
this command. This is because the data is scaled, thinned 
for economy, checked for rotational direction, and clipped 
ao iv els (read) ih: 

7. THINning allows the user to specify the three node 
elimination parameters. These are minimum segment angle, 
minimum segment length, and maximum segment length. All 
three are defaulted in such a way that the data is not 


modified. 
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C. LEVEL TWO COMMANDS 


1. AUTOmatic places the program in an “automatic" mode 
of operation. This means that loop connectivity, loop 
linkup, and triangulation will all be attempted 


automatically. 
2. CAP allows the user to generate caps on contour 
loops. 
3. EDIT allows the user to edit loops by using the 
following commands. 
a. ADD node allows the user to add nodes between 
any two adjacent nodes of the existing contour data set. 
b. DELEte node allows the user to delete any node 
Im the contour data set. 
c. LOOP allows the user to select a new loop to 
edit. 
d. MOVE allows the user to change the coordinates 
of an already existing node. 
ae iINSPect allows the user to inspect the contour data 
ome level at a time. 
9. MANUal allows user assistance in the interpretation 
Of ambiguous contour data by using the following commands. 
a. AUTOmatic causes the program to proceed _ to 
automatically triangulate a group of loops. 
b. BRiDge enables specification of bridges between 
two or more loops on the same level which are branching out 


Of a loop (Cor loops) on the other level. 


go 


c. CLEAr causes the program to delete abli=petyea 


that have been generated between the current pair of contour 


levels. Control is then returned to level three. 

d. GUIDe is a last resort option for Uses 
unreasonable data. The idea is to guide the program in 
triangulation by specifying limits between which to 
triangulate. The following commands are used while in 
GUIDe. 

(1) AUTOmatic invites the program to 


triangulate the loops automatically. 

(2) DENSity allows the user to specify the node 
number label density on the display. 

(3) ERASe allows the user to start Over Onm@ene 
current set of loops. 

(4) RENUmber allows the user to redefine the 
location of the number one node on the level "A" sloop] ge. 
this command if the automatic selection of the node number 
one location is unsatisfactory. 

(5) TRIAngulate allows the user to specify the 
limits on loop numbers, at the two levels, between which 
triangulation is to take place. 

e. INCLude allows’ the user to override the 
alogorithm which selects which loop(s) on level "A= comnege: 


to which loop(s) on level "B". 
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f. NEXT indicates that the user is satisfied with 
the triangulation process for the current pair of contour 
levels and wishes to proceed to the next level. 

6. POST is identical to AUTOmatic, with the exception 
that as each loop is triangulated it is displayed for 


approval. 
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APPENDIX B 


SOURCE CODE SETStiNnes 


The source code for all programs used» 1nN COngunet iam 
with MOVIE.BYU is listed in this appendix. The source code 
for the MOVIE.BYU programs may be obtained from Brigham 
Youne University, Provo, Utah 

Source code is listed in the following order: 
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